Short-time vs. long-time dynamics of entanglement in quantum lattice models 
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We study the short-time evolution of the bipartite entanglement in quantum lattice systems with 
local interactions in terms of the purity of the reduced density matrix. A lower bound for the purity 
is derived in terms of the eigenvalue spread of the interaction Hamiltonian between the partitions. 
Starting from an initially separable state the purity decreases as 1 — (t/r) 2 , i.e. quadratically in 
time, with a characteristic time scale r that is inversely proportional to the boundary size of the 
subsystem, i.e., as an area-law. For larger times an exponential lower bound is derived corresponding 
to the well-known linear-in-time bound of the entanglement entropy. The validity of the derived 
lower bound is illustrated by comparison to the exact dynamics of a ID spin lattice system as well 
as a pair of coupled spin ladders obtained from numerical simulations. 

PACS numbers: 03.65.Ud, 03.67.Mn, 05.50.+q, 02.70.-c 



Introduction 



Motivated by the question whether the time evolution 
of interacting quantum systems can be efficiently sim- 
ulated with the help of matrix-product decompositions 
of the many-body wave function [lH3J| or corresponding 
analogues in higher dimensions [H, HHsl ; the dynamics of 
entanglement in quantum lattice models has become an 
important research area in quantum physics. A conve- 
nient measure of entanglement are @ the von Neumann 
and Renyi entropies. It was shown in [l(| and [ll| that 
the von Neumann entropy of a subsystem which starts 
in an initially separable state has an upper bound that 
grows linear in time. The linear growth of the entropy, as 
observed by Calabrese and co-workers [lU, EH s which cor- 
responds to an exponential growth of the effective bond 
dimension of matrix-product states (MPS) represents a 
severe limitation for the simulability of the unitary time 
evolution of quantum many-body systems. In the present 
note we derive an upper bound to the bipartite entan- 
glement that also holds for short times. In particular we 
consider the purity of the reduced density matrix of one of 
the partitions. General quantum mechanical arguments 
suggest that the purity cannot decease exponentially for 
short times as implied by the linear entanglement bound 
[Tol . ITT1 . [l4| , but rather quadratically. Here we derive a 
quadratic lower bound for the purity. This finding which 
is the main result of our paper has practical relevance 
for the numerical simulation of another class of dynami- 
cal problems that gained a lot of interest recently where 
the time evolution is non-unitary due to a coupling to 
external reservoirs (l5l - [25j . The non- unitary Liouvillian 
dynamics of the system density matrix is equivalent to 
a time evolution of the many-body wave function with 
a complex Hamiltonian and with a stochastic sequence 
of projections, called quantum jumps (26H28j . If the fre- 



quency of such projections is sufficiently large they can 
prevent the growth of entanglement within the system by 
a mechanism similar to the well-known quantum Zeno ef- 
fect [29|. As a result of this the time evolution of open 
system may be simulated using an adaptive MPS expan- 
sion as e.g. within the time evolving block decimation 
algorithm (TEBD) @, i] for longer times. The critical 
frequency of such a Zeno effect for entanglement is de- 
termined by the coefficient of the quadratic term in the 
short time expansion of the purity. Note that this ef- 
fect would be absent for an exponential time dependence 
of the purity. For larger times we derive an exponential 
lower bound corresponding to the well-known linear-in- 
time bound of the entanglement entropy pjj [HI [H| • Al- 
though this result follows directly from the latter entropy 
bound, we derive it here in a few lines without making 
use of the somewhat more involved proof of the linear 
entropy scaling. To illustrate the validity of our findings 
we discuss as an example the ID spin-i XX and XXZ 
models, where we calculate the time evolution of the pu- 
rity using the numerical time-evolving block decimation 
algorithm 2, 3], as well as two coupled spin chains, which 
allow us to illustrate the scaling with the boundary size. 



II — Short-time behavior 

We here consider a lattice model with a bipartition 
into parts A and B, which are, say, compact sets of lat- 
tice sites. We assume that the Hamiltonian of the total 
systems can be written as 



Hab — Ha + Hb 



q e{A\B} 



H^®H^\ 



(1) 
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where Ha and Hb are the parts of the Hamiltonian act- 
ing on lattice sites inside the respective partitions. The 
last sum describes the interaction between the two parts 
an extends over all bonds, labeled by the index q, that 
connect sites from both partitions. We assume an inter- 
action that has strict finite range. In this case the total 



2 



number of bonds scales with the size N of the surface 
separating the two partitions. 

Any pure state of the total system can be decomposed 

as 



fies the equation of motion 



dt 



-itr B < 



Hab , Pab 



(3) 



I*(*))ab = £>/£ 



(2) 



where 



are orthonormal sets of states of 



(4 A) \ and 

the subsystems, £ Q > are the Schmidt coefficients, and 
L is at most the dimension of the Hilbert space of the 
smaller subsystem. As |$ (t)} AB is normalized, = 

1. The reduced density operator of the subsystem A, 
p A = ti b {pab}, where p A B = \^ {t)) A B (* (*)Iab> satis ~ 



(Note that we have set h = 1 throughout this text.) The 
form of the Hamiltonian (fTJ) and the cyclic property of 
the trace allows us to obtain the following equation 



— trp A = -2itr < p A tx B 



(4) 



for the purity rate. The traces on the right hand side can 
be calculated in the eigenbasis of the pa and p B '■ 



tr pAtv B J2 H ( q A) ® ■ PAB )} = £ & tr ^^ A) ® #F W 



/.(A) 



q.cx 



(5) 



q,a,a' 



B) 



(6) 



and in the same way, the second term of the commutator is 
tr <Ltr B Lab ■ E H { q A) ® £< B) J } = E ^ 

I \ 9 /J <j,a,a" 



q 



r 



(7) 



Combining eqs. (|4]), (|6]) and ([7]) one obtains the following and 
differential equation for the purity 



d 
dt 



trp 2 A = -2* e ve^ (&> - oo 



(8) 



ff(B) 



This equation can be rewritten in a compact form 
d 



where 



-^trp^ = tr [6 • Q] , 



(9) 



(10) 



/a/3 



E 



H, 



(A) 



H, 



(B) 
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The matrix Q has only two non-zero eigenvalues. Indeed, 
can be written as 



where 



9 = -2i\a) (b\+2i\b) (a\ , 



\*) = (ti / \ti / \...e I ! 2 ) 
\b) = (e l '\e 2 , \...e I ! 2 



(12) 
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It is easy to show that the nonzero eigenvalues of 6 are 
A± (9) = ±2y/(a \a) (b\ b) - \(a\ b)\ 2 (13) 



Let \q±) be the corresponding eigenvectors of 0, then the 
trace © can be evaluated in the eigenbasis of which 
yields the following equation for the purity rate 



dt 



trp\ = 2yJtipP A -(tT P 2 A f((q + \ Q \q+) - <«-| Q |<?-}) ■ 

(14) 

The right side of this equation can be bounded from 
above by the spread of the eigenvalues of the interaction 
Hamiltonian between partitions A and B: 



Atrpi <2yjtrp\- (tr p\f 



(15) 



due to the second factor the dynamics of the purity does 
not only depend on the local Hamiltonian, but also on 
the purity of the initial state (see examples). 

We can divide the left- and right-hand side by the 
square root term, which after integration yields 



Using 



,,1/4 

— arcsin [ (trp A 



2\l/4 



_d_ 

d7 



(tr^(0)) V4 



(19) 



dr 



(20) 



gives a solution of cq. ([TBI) with the initial purity trp A (0) 



eigenvalues of Y, q e{A\B} ® H^'. 
In a similar way one can show that 



f t trp 2 A > -2^trp A - (trp A Y 



(16) 



Combining the inequalities (|15|) and (jT6j) we obtain 



< 2psjtrp\ - (trp A )\ (17) 



where 



P = K 



\ q 



is a constant that scales linear with the size N of the 
surface separating the subsystems. 

In order to solve the differential inequality (|T7| we 
need an expression or at least an estimate for trp^ in 
terms of the purity. With the help of Hardy's inequality 
(Efc a k) 1,m < (Efe O 17 ' 1 for any a k > 0, and m > n > 
0, (see one can show that trp^ < (tr/r^) 3 ^ 2 . With 

this we find the following differential inequality for the 
purity 



df 



trp A 



< 2/iVWi) 372 - (trp 2 A f- 



(18) 



In general the short time behavior is local, and eq. (JT8J) 
partially supports this intuition, since the control param- 
eter is p, the spread of the local Hamiltonian. However, 



(~|* + arcsin (tr^ (0)) 1/4 , o) 



<tvp A 



< sin 



mm — t 

2 



+ arcsin (trp^(0)) 1/4 , 1)] . (21) 



If trp A (0) = 1, i.e. if the subsystems are uncorrelated 
in the beginning, the upper bound is trivial. The lower 
bound reduces to 



2 4 P^ 

trp A > cos —, 



for pt < 7T. 



(22) 



We note that the lower bound becomes zero at pi > ir 
i.e. it reduces to the trivial one. 

The lower bound can be slightly improved, if we know 
the maximum Schmidt rank L max of the bipartite decom- 
position that can arise along the evolution. One finds 



tvp 2 A > cos 4 lj 



1 
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In 



(23) 



Evidently trp A > L~ n \ x as it should be. Typically 
L max 3> 1 and thus the second term in eq. (j23|) is small. 
In the general be the dimension of the 

smaller Hilbert space. So one must assume this, if there 
is no better bound known a priori. In certain special 
cases as for the ID quantum Ising model with an initial 
state that factorizes in all sites, L max = 2, and the lower 
bound (f2"3"| becomes exact, see below. 



The estimation (|23| follows from the inequality 



4 



trp 6 < 



LI 



(l + V^max - 1) (imax ' trp 2 - 1)^ + 



imax — 1 — V (imax — 1) (imax ' ifp 2 — 1) 
(imax — 1) 



(24) 



which can be proven by the method of Lagrange multi- 
pliers. 

Eqs. (|2"Tj) and (|2"2"|) . resp. (|2"3"|) . provide an estimate 
for the purity for short times. As expected from gen- 
eral quantum mechanical arguments the lower bound 
of the purity decreases quadratically in time following 
~ — (i/r) 2 . The characteristic time r that defines the 
range of validity of the quadratic time scaling is inverse 
proportional to the eigenvalue spread p and scales lin- 
early with the size of the surface N separating the sub- 
systems and thus has an area-law behavior. In order to 
test the scaling of the lower bound with the boundary 
size of the system we consider next a pair of linear spin 
chains subject to an Ising interaction 

N 

H = 1 £d*®Tf, (25) 

3=1 

where a and f are the Pauli operators of the two chains 
respectively. If the initial state is 

\uu) hu-.-iu® in...t) B (26) 



i.e. where all spins in chain A are in the eigenstate of a z , 
resp. f z , with spin down and all spins in chain B in the 
corresponding spin up state, the purity is given by 



tr/? 2 = (cos 4 1 + sin 4 1) (27) 

This can be seen easily, as for the initial product state 
the total purity can be represented as 



N 

trp 2 = trj]p 2 = (trp 2 ) N = (cos 4 1 + sin 4 1) N (28) 

i=l 



because each spin pair of the coupled chains has a purity 
equal to cos 4 t + sin 4 1. We thus see, that the characteris- 
tic time r scales as V^V for this choice of the initial state. 
However, if one considers a Grcenberger-Horne-Zeilingcr 
type initial state 



|GHZ) A ®|GHZ) B = 



71 



(29) 



with |±) denoting eigenstates of a x , resp. f x , one finds 

trp 2 = cos 4 Nt + sin 4 Nt, (30) 

as can be shown by simple algebraic calculations. We see 
that (f30]) coincides with our estimation (f23| since /i = 2N 
and L max = 2. In other words, our estimate is a tight 
lower bound for the purity. In this case the system size 
scaling of the characteristic time is N. 

Note, that in the special case of an Ising Hamiltonian, 
our two-chains example in fact represents also higher di- 
mensional Ising lattices of arbitrary size. Since all sum- 
mands in the Ising Hamiltonian, say on a hyper-cubic 
lattice commute, they can be absorbed in the quantum 
states of the subsystems A an B. They play no role for 
the entanglement, which is only created by the terms di- 
rectly coupling A and B. But those can always be written 
in the form (|25|) . no matter what the spacial dimension 
of the surface is. 



For short times, the quadratic estimate (f22|) is much 
better than any exponential one (|3"4"|) (sec below). Fur- 
thermore it is of fundamental importance. It shows e.g. 
that a sequence of frequent projections of the system onto 
non-entangled states (i.e. no entanglement within the 
system) at a rate larger than p will prevent the build-up 
of such an entanglement in full analogy to the quantum 
Zeno effect [II]. 

Ill — Long time behavior 

The long time behavior of the purity can be obtained 
from the known upper linear-in-time bound of the en- 
tropy [H [H 0. Indeed, by using the convexity of 
— In a;, we immediately get 

s = -I> ln & ^ - ln E £ = - lntr ^< ( 31 ) 

i i 
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Wc thus have 

t r /°A > cxp(— S) > exp (— cq) exp (— c\i) , 

where S < Co + c%t, cq and ci being positive constants. 

In the following we show that (|31j) can also be obtained 
in a simple way from our approach without the necessity 
to invoke the rather involved proof of the linear-in-time- 
bound of the entropy. 

In order to find a suitable estimate for the long-time 
behavior of the purity one has to find a different way to 
bound the right hand side of eq. (JEJ. The interference 
effects become negligible at t ^> ^ and therefore wc may 
use inequalities of the Schwartz type. In other words, we 
may sum all interactions (matrix elements) in modulus 
instead of amplitude. In this case one finds 



d 2 
dt tTPA 



<2j2 



H(B) 
i 



(32) 



H. 



(A) 



H 



of' 



Here we have used 



Making use of Schwartz's inequality, we obtain 




We thus arrive at 



dt 



trp 2 A 



(33) 



where 



V^2 'y ^ \ A max (^Hq M A max 



The solution of this differential inequality with the initial 
condition trp\ (0) = 1 is 



trp^ > exp(-xi). 



(34) 




FIG. 1: (color online) Time evolution of the purity in the 
80 + 80 site XX chain compared to the bounds (J22j| and (|34)) . 
see text. The inset shows a closeup for short times. Note 
that t has no units since we have choosen the dimensionless 
Hamiltonian (1351) - 



IV — Numerical examples: spin— i XX- and 
XXZ-chains, multidimensional quantum Ising model 

To illustrate the quality of the bounds given in (f2"2"j) and 
P4p we perform exact numerical simulations for a large 
ID spin system as well as for two coupled spin chains of 
small size. We first consider the spin-i XX-model 



(35) 



For this model the constants which enter our estimates 
take on the values fi = 2 and x = V%- We look at a chain 
of 160 spins and in order to maximize the dimension of 
the subsystems we choose an equal partition with A (B) 
being the left (right) half chain. The initial state of the 
system is taken to be a product state, specifically the 
anti-ferromagnetic state 

l* A (f = o)> = i* B (t = o)) = \mm • • • n> ■ (36) 

We choose this particular initial state in order to have a 
large maximum entropy since I^a^ = 0)) corresponds to 
half filling, so the dimension of the Hilbert space accessi- 
ble with respect to the present conservation of the total 
z-magnctization is maximized. The purity is initially 
1. Fig. [TJ shows the evolution of the purity over time. 
The red, straight, thick line is the numerical results from 
our simulation using the time-evolving block decimation 
(TEBD) method @,[|. This results can be considered 
numerically exact as discussed below. The solid, black, 
thinner lines show the bounds (f2"2"j) . The one starting at 
t = indicates that our bound is optimal up to second 
order for times short compared to the inverse coupling 
between A and B, if we start initially from a pure state. 
However when starting from an initially entangled state, 
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FIG. 2: (color online) Time evolution of the purity in the 
80 + 80 site XXZ chain for A = | compared to the bounds 
(|22[) and (|34[) in analogy to Fig. Q] The inset shows a closeup 
for short times. The dotted lines show the purity in the XX- 
model (Fig. [l} for comparison. 

we can only expect to get agreement up to zeroth order 
from (f2Tj) . as illustrated by the black, solid, thinner lines 
starting at t = 1 and t = 2. While the exponential lower 
bound plotted in the dashed line, is bad for short 

times, it has the property of remaining finite for all times 
in contrast to (f2"2"j) . So one can smoothly concatenate the 
two bounds at time 

t\ = — arctan [ — ] , (37) 

assuming we started with a pure state at t = 0. This 
combined bound is superior to both the quadratic short- 
time and exponential long-time estimates and is shown 
as a green, dot-dashed line. 

Analogous calculations where also done for the spin-^ 
XXZ-model 

j 

choosing the anisotropy A to be |. This again yields 
a constant of [i = 2 for the short time behavior. But 
while the exponential bound increases to \ = and 
one could expect a faster decay of the purity due to the 
fact that this system can not be mapped to free fermions, 
the true curves are very much alike for both systems, see 
Fig.H 

The numerical calculation was done using a matrix di- 
mension of D — 500, a time-step width of 0.01 in a fourth 
order Trotter decomposition, and exploiting the conser- 
vation of the total magnetization explicitly. Although the 
dimension of the subsystems is adaptively truncated to 
D, this cannot introduce error on the timescale plotted, 
since the purity remains well above the minimal value of 
l/D representable using matrix product states. Thus all 
relevant states are included by the algorithm. 



FIG. 3: (color online) Short time scaling of the purity with 
boundary size for a pair of coupled spin chains with Ising in- 
teraction. Exact diagonalization was done for t = 0.001, such 
that the quadratic order in (|22[) and (J2SJ suffices to describe 
the result. While the GHZ-like initial states show the fastest 
increase of entanglement (r ~ N, eq. Q22p . straight line), and 
the product like ones stick to the moderate r ~ y/N (eq. (|28[) . 
dashed line), W-type states show an intermediate behavior. 

Figure [3] illustrates the tightness of ([22]) for proper ini- 
tial conditions. It shows 1 — trp^ in a system of two Ising 
spin chains of length N, subject to the Hamiltonian (|25|) . 
after a short time of t = 0.001 of evolution. Ising-typc 
couplings inside the chains where also taken into account, 
but do not contribute to the entanglement between the 
two chains. While for simple product states between the 
sites, one can expect a scaling of r ~ (dashed line), 
we know from (|30p , that there are in fact initial states, 
that contain sufficient entanglement along the boundary, 
to giver ~ N, i. e., where (|22[) is tight. The different 
symbols correspond to different initial states, the evolu- 
tion of which was calculated via an exact diagonalization. 
As already seen above, we have a N scaling for an initial 
product state like (|26p . while we get an exact N 2 scal- 
ing for GHZ-like initial states (|29|) . Also shown are other 
states like the W-type states |Wj,) A ® \ W P ) B , where 

\W P )= Y, & x h al...&l\H...i). (39) 

l<ji<h<-<j P <N 



Summary 

In summary we derived an upper bound for the time 
evolution of the bipartite entanglement in quantum lat- 
tice models in terms of a lower bound to the subsys- 
tem purity. As one would expect from general quantum 
mechanical considerations the purity decreases for short 
times quadratically in time. The corresponding charac- 
teristic time was shown to be limited by the spread of the 
eigenvalues of the part of the Hamiltonian that accounts 
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for the interaction between the partitions. The latter 
scales linear with the size of the surface separating the 
two partitions and thus the entanglement follows an area- 
law behavior. For larger times we derived a lower bound 
of the purity that decrease exponentially in time. The 
latter is equivalent to the known linear increase of the en- 
tanglement entropy in the long-time limit. The existence 
of a quadratic short-time bound means that a sufficiently 
frequent sequence of projections to non-entangled states, 
as for example due to a dissipative process, can prevent 



the build-up of entanglement within the lattice system. 
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